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ANALYSIS OF DEGENERATE CROSS-DIFFUSION POPULATION 
MODELS WITH VOLUME FILLING 

NICOLA ZAMPONI AND ANSGAR JUNGEL 


Abstract. A class of parabolic cross-diffusion systems modeling the interaction of an arbi¬ 
trary number of population species is analyzed in a bounded domain with no-flux boundary 
conditions. The equations are formally derived from a random-walk lattice model in the 
diffusion limit. Compared to previous results in the literature, the novelty is the combina¬ 
tion of general degenerate diffusion and volume-filling effects. Conditions on the nonlinear 
diffusion coefficients are identified, which yield a formal gradient-flow or entropy structure. 
This structure allows for the proof of global-in-time existence of bounded weak solutions 
and the exponential convergence of the solutions to the constant steady state. The existence 
proof is based on an approximation argument, the entropy inequality, and new nonlinear 
Aubin-Lions compactness lemmas. The proof of the large-time behavior employs the en¬ 
tropy estimate and convex Sobolev inequalities. Moreover, under simplifiying assumptions 
on the nonlinearities, the uniqueness of weak solutions is shown by using the H~^ method, 
the E-monotonicity technique of Gajewski, and the subadditivity of the Fisher information. 


1. Introduction 

In this paper, we analyze a class of mnlti-species population cross-diffusion systems with 
volume-filling effects. Such systems arise in various applications, like spatial segregation of 
interacting species [25] . chemotactic cell migration in tissues |24j . and ion transport through 
membranes [8]. Our model class can be derived from a system of random-walk master equa¬ 
tions in the diffusion limit for a large class of transition rates (see Appendix . The key 
novelty of our analysis is the identification of a new entropy or formal gradient-flow structure 
and the treatment of non-standard degeneracies in the diffusion coefficients, which signifi¬ 
cantly extends previous results in m- 
The diffusion systems have the form 

(1) dtu — div(A(u)Vn) =0 in 0, t > 0, 
with boundary and initial conditions 

(2) {A{u)Vu) ■ u = 0 on (90, t > 0, n(0) = in Q. 

Here, II C (d > 1) is a bounded domain, A{u) = {Aij{u)) £ is a diffusion matrix, 

the function u = {ui,... ,Un) '■ H x (0,oo) —>• M"" is the vector of the proportions of the 
subpopulations, and Un+i = 1 — proportion of unoccupied space. In particular. 
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0 < < 1 for all i = 1,... ,n + 1. The ith component of equations ([T]) and m has to be 

understood, respectively, as 


n n 

dtUi — ''^^div{Aij{u)Vuj) = 0 , Aij{u)Vuj ■ v = 0. 

i=i i=i 


The boundary condition in ([2]) means that the physical or biological system is isolated; the 
species cannot move through the boundary. For ease of presentation, we have neglected 
reaction and drift terms in the equations. We refer to Section [3 for a discussion of more 
general models. 

The diffusion matrix in (JT]) is given by 


dpi 


(u), i,j = l,. 


, n, 


(3) Aij{u) = 6ijPi{u)qi{un+i) + UiPi{u)qi{un+i) + Uiqi{un+i) 

where 6ij is the Kronecker delta. The nonnegative functions pi and qi model the transition 
rates in the random-walk lattice model. The coefficients Aij are derived from this model in 
the diffusion limit (see Section The function q^ vanishes when the cells are fully packed, 

i.e. if X]r=i ~ 9i(0) = 0 a-iid qi is nondecreasing. In the literature, several special 

models were considered and we review now some of them. 


Example 1. 1. Population-dynamics models. The case n = 2, Pi{u) = Oio + anui + 0 , 2^2 
and qi{u 3 ) = 1 for i = 1,2 was suggested by Shigesada, Kawasaki, and Teramoto [25] to 
describe the spatial segregation of interacting populations and to study the coexistence of 
two similar species. This model has attracted a lot of attention in the literature. One of the 
first existence results is due to Kim |2r)] who imposed some restrictions of the parameters 
aij. The tridiagonal case 021 = 0 was investigated, e.g., by Amann [T] and Le [2T]. The first 
global existence result without any restriction on the diffusion coefficients (except positivity) 
was achieved in [T8| in one space dimension and in [U |T0| in several space dimensions. The 
case of concave functions pi and p 2 was analyzed by Desvillettes et al. m, recently improved 
in m- The n-species case with superlinear functions pi (u) was investigated in [T9] ; also see 
[3] for a so-called relaxed system. 

2. lon-transport models. The case Pi{u) = 1 for z = and q{un+i) = Un+i was 

employed to describe the motility of biological cells m or the ion transport through nanopores 
[8]. The global existence of bounded weak solutions was proved in [?]■ This result was 
generalized in |19j to a class of nondecreasing functions including all power functions q{s) = 
with a > 1. The models in [HJIJT] also include a drift term to account for electric effects, and 
we discuss these extensions in Section [3 

3. Multi-species chemotaxis models. A special case of the model in [ 23 ] is given by Pi{u) = 1 

and q{un+i) = itn+i; similar to the ion-transport model. In fact, the system in |24j contains 
additional terms which cannot be described by ([3]) since the transition rates assumed in |24j 
are not of the type pi{u)qi{un+i) (see (l66l) in Appendix lAl) but they equal Pi{u) + qi{un+i). 
We refer to the discussion in Section [3 D 


In the model classes (i) and (ii), either pi = 1 or q^ = 1. In contrast, we investigate here 
a more general model class allowing for nonconstant functions pi and q^. A guiding example 
is system ([T]) with diffusion coefficients ([3]) and Pi{u) = uiU2, qiis) = s for i = 1,2, which 
models volume-filling effects in population systems. The diffusion matrix reads explicitly as 


(4) A{u) = 


rii(l - ui - U 2 ) + (ui -I- U2){1 - U 2 ) 
U2 


Ul 


U2{1 - Ul - U 2 ) + {ui U2){1 - Ul) 
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We will show in Theorem [T] that ([T]) with this diffusion matrix possesses a global weak solution 
satisfying 0 < u{t) < 1 for all t > 0. In fact, Theorem [T] is concerned with much more general 
models. 

The analysis of system ([T|) with diffusion matrix ([3]) faces a number of mathematical chal¬ 
lenges. First, the equations are strongly coupled such that standard tools, like maximum 
principles and regularity theory, generally do not apply. Second, the diffusion matrix is gen¬ 
erally not positive definite and thus, even the local-in-time existence of solutions is nontrivial. 
Third, since the variables Ui are proportions, we need to prove lower and upper bounds for the 
solutions (here, Ui >0 and Ui < 1), but maximum principle or invariant region methods 
seemingly do not apply. Fourth, the parabolic system may be degenerate (e.g. like in Q for 
u = (0,1) or ri = (1,0)). 

Some of these difficulties have been dealt with in, e.g., m under the assumption that the 
diffusion system has a formal entropy or gradient-flow structure, i.e., there exists a convex 
functional h : IF —)■ (called entropy density), where ^ C M”, such that the matrix B = 
A{u)h''{u)~^ is positive semi-definite and ([T]) can be written as 

(5) dtu — div(i?V/i'(tt)) = 0, 

where h'{u) and h”{u) are the Jacobian and Hessian of h, respectively. This formulation has 
two advantages: First, H[u] = h{u)dx is a Lyapunov functional along solutions u{t) to 

O-®, 

(6) ~^[^(^)] = [ h'(u) ■ dtudx = — f Vu : h"(u)A{u)'Vu = — f Vw : BVwdx < 0. 

dt Jq Jq Jq 

where w = h'(u) are called entropy variables. In particular, this yields a gradient-type 
estimate for w or u. Second, if h' is invertible on (see Lemma [5|), the original variable 
u = {h')~^{w) is an element of Thus, if ^ is a bounded domain, we obtain lower and 
upper bounds for u without the use of a maximum principle. In our situation, we define 
^ = {u G M” : Uj > 0 for f = 1, ..., n, such that Ui is positive and bounded by 

one. 

There remain still two issues for systems with diffusion coefficients ([3]). The first one is 
to identify a suitable entropy density h, the second one is the possible degeneracy. In the 
example given by ([1]), we choose 

2 

h{u) = '^ufilogui - 1) (1 - ui - rt 2 )(log(l -ui- U 2 ) - 1) 

i=l 

+ (ui + U 2 )(log(ui -I- U 2 ) - 1) -I- 4, 


which yields the matrix 
B = A{u)h''{u)-^ 


Ul{ui + U2){1 - Ml - M 2 ) 
0 


0 

M2(mi -I- M 2)(1 - Ml - M2) 


At least one eigenvalue of B vanishes if m G dS> = {mi = 0, M2 = 0, 1 — mi — M2 = 0}. In 
this sense, system m is called to be of degenerate type. Generally, systems m are always of 
degenerate type since g(0) = 0. Here, we develop a technique to deal with such a degeneracy. 

We overcome these issues by developing two main ideas. Our first key idea is the identih- 
cation of a class of functions pi and qi for which we are able to dehne a novel entropy density. 
The second idea is the extension of the Aubin-Lions compactness lemma to non-standard 
degenerate cases. In the following, we detail these concepts. 
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We make the following structural hypotheses on the functions pi and qp There exist func¬ 
tions g : [ 0 ,1] —>■ M, X : ^ M and a number 7 > 0 such that for alH = 1 ,..., n, 

(7) q{s) := qi{s) > 0, g'(s) > 79 ( 5 ) for s G (0,1), g(0) = 0, g G ^^([0,1]), 

( 8 ) Pi{u) = exp ( —- ) for M G IF, X > is convex on x ^ C^(^). 

\ J 

Examples of functions q and pi satisfying these conditions are given in Remark [2j We define 
the entropy density 


(9) 


^ ^ rUn+l 

h{u) = '^{uilogUi-Ui + l)+ log q{s)ds + xiu), uG 

i=i 


where a G (0,1] is such that jl^logq{s)ds > 0 for all b G (0,1), namely 


( 10 ) 



if g(l) < 1, 
ifg(l) > 1. 


Notice that we require that all functions qi are the same and that pi possesses a particu¬ 
lar structure. It seems to be difficult to treat more general cases, except imposing other 
conditions. 

Surprisingly, system ([I]) with ([3]) partially decouples in the entropy variables. Indeed, we 
may write the following formal “generalized” gradient-flow formulation 

( dh \ 
g('Un+i)^Vexp — j = 0 , i = l,...,n, 

which makes the degenerate structure more apparent than ([I]). We also note that \i q = 1, 
we obtain dtUi = A(exp(9/i/9Mj)) = A{uiPi{u)). This structure was exploited in [THlfT^ . 

A computation, which is made rigorous below, shows that the following entropy inequality 
holds: 


where c > 0 is some constant. We wish to deduce gradient estimates for tti,..., which 
are needed to apply the Aubin-Lions compactness lemma for a suitable approximated system. 
However, because of the degeneracy of q (i.e. q{0) = 0), these estimates are nontrivial. We 
overcome this problem by proving two compactness results. 

The first compactness result essentially states that if we have (i) uniform gradient estimates 
for the bounded sequences (^g) and (^e%), (ii) a uniform estimate for the (discrete) time 
derivative of and (iii) the strong convergence ^ in L^, then up to a subsequence, 

^efiVe) Cfiv) ill for any continuous function / (Lemma[7|). If were strictly positive, 
the statement would be a consequence of the usual Aubin-Lions lemma [26]. Here, we are 
able to deal with functions which may vanish locally. The case f{s) = s was considered in 

|71[I9|. 

The second compactness result is a generalization of the Aubin-Lions-Dubinskh lemma; 
see, e.g., [IT1[22]. It states that if a bounded sequence (ue) possesses a uniform estimate for 
the (discrete) time derivative and a uniform gradient estimate for ( 5 (tie) and Q'{us) for some 
nonnegative convex increasing function Q, then up to a subsequence, ^ u strongly in 



q{un+i)'^ ^ \Vq{un+i)^^'^\'^ I dx < 0, 


2=1 
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(Lemma [ 8 |). This result is complementary to the nonlinear Aubin-Lions lemma stated in [ 22 ] 
and generalizes the lemma in m stated for Q{s) = s" with a > 1. 

Based on the above ideas, we prove three results. First, we show the global-in-time existence 
of bounded weak solutions to satisfying the entropy inequality ([IT]) (Theorem [T|). 

Second, the entropy inequality and a convex Sobolev inequality allow us to show that Un+i{t) 
converges to the constant steady state in the sense. Moreover, if q is strictly positive, this 
convergence also holds for ui{t),, Un(t) (Theorem [3|) . Third, if pj = 1 for alH = 1,..., n, 
there is a unique weak solution to ([I])-(l3|). The proof combines the H~^ method and the 
Fi-monotonicity technique of Gajewski m- 

The paper is organized as follows. The main results are stated and commented in Section [2j 
Section [3| is devoted to the proof of some auxiliary results, like the positive semi-definiteness of 
the matrix h''{u)A{u) and the Aubin-Lions compactness lemmas. The three main theorems 
are proved in Sections 0] (5] and (U respectively. Extensions of our model are discussed in 
Section [T] Appendix A is concerned with the formal derivation of ([T|) from a random-walk 
lattice model. 


2. Main results 

We state our main theorems and detail the ideas of the proofs. The hrst theorem is 
concerned with the global existence of bounded weak solutions. Recall that 

n 

(12) ^ = (ui,..., Un) G IR"" : Uj > 0 for i = 1,... ,n, Uj < l|. 

i=i 

Theorem 1 (Global existence). Let T > 0, let : Ll ^ ^ be a measurable function such 
that h{u^) G and let A{u) he given by ([3|). Assume that hypotheses ([7|) and dH) hold. 

Then: 

(i) There exists a weak solution n : x (0,T) —^ to ([Il)-(l2|) satisfying Ui > 0, Un+i ■= 
1 - Ya=i ^ 0; 

(13) g(u„+i)^/^, UiPi{u)q{un+if^‘^^L^{{),T-,H^{Ll)), 

(14) Ui G L“(0,T;L°°(O)), dm e L^{{),T]H^{n)'), f = l,...,n. 

The function u satisfies the weak formulation 

" T n j-T j- 

{dtUi,(l)i)dt + '^ / / [q{un+i)^^‘^V{uiPi{u)q{un+if^‘^) 

Jo Jn 

- 3uiPi{u)q{un+iy^'^"^q{un+i)^^^] ■ Vcfidxdt = 0 

for all (fi,... ,4>n G L^{0,T; and u(0) = in the sense of . Here, (•, •) 

denotes the duality product of H^{Qy and 

(ii) The following entropy inequality holds: 

(16) [ h{u{t))dx + co [ [ (^q{un+i)‘^\Vuy‘^\‘^+ \Vq{un+i)^^‘^\A dxdt 

Jn Jo Jn J 

< [ h{u^)dx, 

Jn 

where cq = 4po min{l, (5} > 0 with po and 6 being defined in ([25|1 below. 
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(iii) If /q I log <7(5) Ids = +00 for all 0 < b < 1 then tt^+i >0 a.e. in Q x (0,T). 

Remark 2. We present examples of functions q and pi satisfying d?]) and ([ 8 ]), respectively. 
Hypothesis d?]) is satisfied by q{s) = s" for s G [0,1], where a > 1. Indeed, the inequality 
Q'is) > jqis) holds for all s G [0,1] with 7 := a. Another example class is given by q{s) = 
exp(/(s)) — 1 with /(O) = 0 and f'{s) > 7 > 0 for s G [0,1]. A concrete example is 
q{s) = exp(s“) — 1 with 0 < a < 1. A third example is q{s) = exp(—with a > 0 which 
satisfies the assumption stated in Theorem [H part (iii). 

Hypothesis dS]) is satisfied by every function pi{u) = pi{ui), where pi G C'^([0, 1 ]) is strictly 
positive and nondecreasing. Indeed, let us define 



for s G [0,1], i = 1,..., n, and u = (ui,..., Un) G Here, k > 0 is such that Xi > 0 in [0,1]. 
Since pi is strictly positive and nondecreasing in [0,1], it follows that x"(u), given by 


d'^X 

duiduj 


(u) 



i,j = 1 ,... ,n, 


is positive semi-dehnite and x'- ^ ^ [0,oo) is convex. Furthermore, exp(clx/fluj) = Pi{ui) = 
Pi{u) for u G 

Another example is given by Pi{u) = with a* > 0, i = 1, ..., n. Indeed, the 

function x('w) = convex on ^ and satisfies exp{dx/dui) = 

exp{ai ajUj)) = Pi{u). This example corresponds to the diffusion matrix (jH) for 

n = 2 and oi = 02 = 1. □ 


The proof of Theorem [T] is based on an approximation and regularization of dH) • More 
precisely, we consider the semi-discrete system 

— {u{w'^) — u{w^~^)) = div{B{w'^)Vw’^) + + w'^) + £t‘^ (— 

r ^^ 

2<\a\<m 


with homogeneous Neumann boundary conditions, where t > 0, £ > 0, m > d/2, u{w) = 
{h')~^{w), approximates w{kT), and is a partial derivative of order 2|a|, with a G Nq 
being a multiindex. Compared to |19] . we need two regularization levels: the regulariza¬ 
tion given by Aw^ -|- and the regularization given by the sum over a. The second 
regularization is needed to obtain approximate L°° solutions (observe that 
while the hrst one allows us to interpret the weak formulation in the larger space instead 
of This is needed to apply the generalized Aubin-Lions Lemmas [7] and [H for which 

H~^ is required. 

The entropy inequality (fTTI) . adapted to the above problem, yields uniform estimates. 
Hence, applying the Leray-Schauder fixed-point theorem, we obtain the existence of semi¬ 
discrete solutions. The same entropy inequality provides a priori estimates uniform in r 
and £. First, we perform the limit e ^ 0, then the limit r —)• 0. The latter limit is highly 
nontrivial since we have only an bound for q{un+i)'^uy‘^, and q{un+i) = 0 at Un+i = 0 is 
possible. This degeneracy will be overcome by the compactness result in Lemma [Tj 
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The second result is about the large-time behavior of the solutions to the constant steady 
state given by 

We are able to prove exponential convergence of Un+i{t) and, under an additional assumption 
on q, also oi ui{t),... ,Un{t). 

Theorem 3 (Convergence to steady state). Let O be eonvex, G let A(u) be 

given by (l3|), and assume that dZD and ([8|) hold. Furthermore, let q G 1]) be such that 

q' is strictly positive and q/q' is concave on (0,1). Let u : Q x {0,T) ^ ^ be a weak solution 
to ([I])-© in the sense of Theorem{^ Then 

( 17 ) \\un+i{t) - u^+,\\L 2 (^n) < Cie-^^\ t>0, 

where Ci = ( 2 / 7 )^/^ ||/i* and Ai = cogi/( 4 c 5 ), h* is the relative entropy density 

(see (fTOll j. qi ;= min^g^o,!]> 0 , cq > 0 is defined in TheoremUi and cs > 0 is the 
constant of the convex Sobolev inequality in Lemma[TU Moreover, if qo := min^gjo^j q{s) > 0, 

( 18 ) \\ui{t)-uf\\L 2 ^n) <Cie~^^\ t>0, i = l,...,n, 

where A 2 = c^qo/cL and > 0 is the constant in the logarithmic Sobolev inequality (see, e.g., 
m Lemma 1\). 


The convexity of Ll and the concavity of q/q' is needed to apply the convex Sobolev in¬ 
equality (see Lemma [To] below). For instance, q/q' is concave for q{s) = s“ with a > 0. The 
condition on the strict positivity of q contradicts the assumption g(0) = 0 in Hypothesis Q. 
However, Theorem [1] is also valid for functions q(0) > 0. In fact, the existence analysis is 
much easier in this case since the problem becomes nondegenerate. 

The idea of the proof is to derive an inequality for the relative entropy 

(19) f h*{u\u°°)dx= f [h{u) — h{u°°) — h'{u°°) ■ {u — u°°))dx. 

Jo. Jo. 

A computation, which is made rigorous in Section (5] shows that 


fL 

dt 



^n + l Cd 


In 


logq{s)dsdx + c / \'Vq{un+i)^''^\‘^dx < 0 

Jn 


for some c > 0. The entropy dissipation can be bounded from below (up to a factor) by the 
relative entropy by means of the convex Sobolev inequality [2]. Together with the Gronwall 
lemma and the convexity of the relative entropy, this yields exponential convergence of Un+i{t) 
to u//(_i in the L^ norm. In a similar way, we obtain the entropy inequality 


/L 

dt 


In 


n 

E 

2=1 


Ui{t)\og 


^^^^dx + c [ q{un+if\Vuy^fdx < 0 . 
Jn 


Here, the degeneracy of q at Un+i = 0 prevents the application of the logarithmic Sobolev 
inequality. For this reason, we assume that q is strictly positive. Then, by Gronwall’s lemma 
again, we deduce the exponential convergence of ufit) to u/° in the L^ norm. 

Our last theorem is a uniqueness result in the special case pi = 1. This includes the 
ion-transport model [ 8 ]. 
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Theorem 4 (Uniqueness of solutions). Let the assumptions of Theorem{l\ hold and let pi = 1 
for i = 1,... ,n. Then there exists a unique weak solution to dU)-® satisfying dH-dll. 

The idea of the proof is to combine the H~^ method and the i?-monotonicity technique of 
Gajewski m- In fact, we exploit the special structure of ([I|) and ([3]) in the case pi = 1: 

dtUi = div (g(un+i)Vui - ufS/q{un+i )), f = 1,..., n. 

Summing all these equations, we end up with a simple equation for Un+i'- 

dtUn+i = AQ(u„+i), Q'{s) = q{s) + (1 - s)(f{s). 

The uniqueness for Un+i is shown by the usual H~^ method. The uniqueness for the remaining 
components Ui is more difficult since we cannot easily treat the drift term. This is in contrast 
to the drift-diffusion equations for semiconductors, where a monotonicity property of the drift 
term can be exploited. Here, we employ the U-monotonicity method m- This method is 
based on the convexity of the logarithmic entropy. More precisely, define the distance 

d{u, v) = + ?(u) - 2^ ^ dx, 

^(s) = s(logs — 1) + 1, s > 0. 

A formal computation, which is made rigorous in Section [6l using the subadditivity of the 
Fisher information (see Lemma [9|), shows that 

■^d(u, u) < 0, t > 0, 
dt 

and consequently, d{u{t),v{t)) < d{u{0),v{0)) = 0 for f > 0. Since ^ is convex, we infer that 
d{u{t),v{t)) > 0, which finally yields Uj = u, for z = 1,... , n. 

3. Auxiliary results 

3.1. Invertibility of the entropy transformation. We show that the transformation of 
variables w = h'{u) can be inverted. Recall that the set ^ is defined in (1121) . 

Lemma 5. Let assumptions (I7|)-(l8|) hold. Then the function h : ^ M, defined in ([9]), is 
strictly convex, nonnegative, belongs to C'^(^), and its gradient h' ^ ^ is invertible. 
Moreover, the inverse of the Hessian h" ■. Q ^ M” is uniformly bounded. 

Proof. We first show that /i' : ^ M” is invertible. For this, we observe that 
— = logUi - logqi^l i = l,...,n. 

* j=i * 

The Jacobian of the function g = {gi,..., g^) ^ ^ M", defined by gi{u) = logUj — logg'(l — 
YlJ=i is positive definite since 

% _ ^ q'jun+l) 
duj Ui q{un+i)' 

It is shown in Step 1 of the proof of Theorem 6 in m that 5 : ^ M"" is invertible. Thus, 

we can define the function f = h' o g~^ : ^ M”. Since h''{u) and g'{u) are nonsingular 

matrices for u G S’, the Jacobian of /, 

ny) = h''{g-\y)){gr\g-Hy)), 
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is nonsingular for y £ M”. Moreover, by the definitions of / and g, we have 

( 20 ) f{y) = y + x{g~Hy)), 2 /eM". 


Hypothesis ([8]) states that x' £ C°(^) C thus (I2U1) implies that \f{y)\ —)• oo as 

\y\ —)■ oo. This property as well as the invertibility of the matrix f'{u) allow us to apply 
Hadamard’s global inverse theorem, showing that / : M” —)• M” is invertible. Consequently, 
also h' = f o g ■. 3i ^ M"' is invertible. 

It remains to prove that the inverse of the Hessian of h is bounded. Since q'/q > 0, 
0 < < 1, and X is convex in the expression 


( 21 ) 


d'^h ^ ^ q'iun+i) d'^X 

duiduj Ui q{un+i) duiduj ’ 


shows that v~^h"{u)v > \v\'^ for all u G u G M"". We infer that all points in the spectrum 
of h" are strictly positive in S>. In particular, h is strictly convex. As h" is symmetric, we 
conclude that the inverse of h” is bounded in 3!. □ 

3.2. Positive definiteness of HA. We show that the product HA of the Hessian H := h''{u) 
and the diffusion matrix A = A{u) is positive definite. This result is needed to deduce gradient 
estimates for tt; see ([6|). 


Lemma 6. Let assumptions ©-([H]) hold. Then the matrix HA is symmetric and positive 
definite. More precisely, for all u € 3 and v G M”’, we have 


( 22 ) 

where 


v^{HA)v > poq{un+i) 'y' — +Po6 

^ Ui 


Q'juur^+i? 

q{Un+l) 



(23) 


Po = min inf Pi{u) >0, 6 = min 

l< 2 <n 


1 , I > 0 

2 snpiy 2 <,<i 9 '(s) J 


Proof. First, we verify the symmetry of HA. Using (1211) and the definition of A, we find that 




k=l 


Ui duiduk q 


dkjPkQ + UkPkq' + Ukq 


dpk \ 
duj ) 


^ pm , dpi d^x d^X 

= —+piq+—q+ g^Pjq + H-p^^kPkL 

LL'!^ L/Lt'j KyLLjjiyLLj 


k=l 

U\2 


duiduk 


d X dpk , {q'Y dpk 


q 


k 


k=l ^ ^ k=l 

Dividing this equation by q, defining (p = q'/q, and taking into account that, by assumption 


(9^X 1 dpj 1 dpi 

duiduj pj dui Pi duj 


for z,j = 1 ,... ,n. 


-{HAfi,-S,,-+Pi^ + — 


+ 


dui 


n rA 

. dpk 


k=l 


we infer that 
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71 (-\ r\ 71 71 r\ 

^ S^ OpkOpkUk ^ , 2 opk 

J 4—1 k=l 


fc=l 


^ Pi ^ dpi ^ dpj ^ ^ 
Ui dui dw ^ 


k=l 


k=l 

Uk dpk dpk 
Pk dui duj 


(24) 


+ ^ ipi+Pj + '^Uk 

V k=l 


dpk dj^ 

dui duj 


+ ^'^'^PkUk, 
k=l 


which proves the symmetry of HA. 

Next, we show the lower bound (l 2 ^ . Since pi is strictly positive in pi{u) = A + pi(u) 
for any A G (0,po)) where po > 0 is defined in (1231) . and Pi{u) is still strictly positive in 
Then we can write (p41) as HA/q = M + \N for two matrices M = (Mjj) and N = (Nij), 
defined by 


., _ . K ^ I ^ I Ufc dpk dpk 

Ui duj dui ^ Pfc + A dui duj 

k=l 


+ <p ipi+Pj + '^Uk 

V k=l 


dpk ^ 

dui duj 




Pk^k •} 


k=l 


Nij = -^ + 2ip + if^'^Uk = — + 2(p + ip^{l - Un+l). 

Ui Ui 


Let V G R"'. Then v (HA/q)v = v Mv + v Nv. We consider v Nv first: 

/ \ 2 


(25) 


The inequalities 


i=l 


j'^Nv = '^'^ + ip{2 + (p{l- Un+l)) ' 


^ 2 = 1 


2q{s) + (1 - s)g'(s) > (1 - s)q'{s) > -^q'{s) 

2q{s) + (1 - s)q'{s) > 2q{s) > - 

supi/ 2 <,T<i q'[(T) 

imply that 

2g(u„+i) + (1 - tt„+i)g'(u„+i) > 5q'{un+i), 

where <5 > 0 is defined in (1231) . Thus, (1^ yields 

no / r) \ ^ 

V 


for 0 < s < -, 


for - < s < 1 , 
2 - - > 




2 = 1 


^ 2 = 1 


Finally, we show that Mv > 0, which, together with the above estimate proves the 
lemma. Using the definition of M, we compute 


(26) 


,,T 


71 ^ 71 

Pi 2 i ^k 




dpj 


2=1 


k=l 


Pk 


^ 2=1 


dui 


*j'=i 




2 
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n n 


I j=l / \i=l k=l i=l 

Let US consider the terms proportional to and 99^ 




\k=l 


vi=i 


n n 




1 ,J = 1 / \i=l 

( n 

'^UkPk 
\k=l / 

/ n ^ 

'^UkPk 
Kk=l j 


' J / ,1 

—^ OUi 

:=1 i=l . 

+2Jt 


\k=l 


u'=i 


<j=l 




J27=iPiVi + ELi Er=i Vijdh/di 

E n ^ 

k=l '^kPk 


P^Vj + 

i=i 


TJi=iPiVi + TJk=i '^k TJi=i Vi{dpk/dui) 


E n 

k=l "^kPk 

Etifi". + SLi Ef-i Vi{Bpk/dui)f 


E n ^ 

k=l '^kPk 

Inserting this expression into (1261) yields 


n / n 


„^Mv>y'!^vf + y'^[y'vM] +2V 


9m 


'3 

a'^i'^3 
OUi 


(Er=iK^*+ELi Er=i 

E n ^ 

k=l "^kPk 

We claim that the right-hand side can be written as a square. To see this, we introduce the 
vectors y = {yi,..., yn), z = ( 21 ,..., z„) G M"- by 


yi = \ —V 


Ui opi 

'i + \ — 

V Pi ^^ 

V /c=l ^ 


9pi _ _ y/ 


Zi = 




\/Efc=l '^kPk 


, z = l,...,n. 


The properties 


n / n 


bi"= 


1 . + +^E 


9pfc 


2 = 1 


Pk dui 


k=l ^ \i=l 


*J=1 


du 


-ViVj, 


Yli=iPiVi + Ylk=i'^kYli=iMdpk/dui) 

y ■ z — 


show that 


■\/Efe=l '^kPk 

^Mv > |yp - {y ■ z)"^ = \y - {y ■ z)z\‘^ > 0 . 


The lemma is proved. 


□ 


3.3. Generalized Aubin lemmas. We prove two generalized Aubin lemmas for functions 
which are piecewise constant in time, extending results from mm- 

Lemma 7 (Generalized Aubin lemma I). Let • • •, iVn^) be sequences of functions 

which are piecewise constant in time with constant step size r > 0 and which are bounded in 
L°°(0, T; L°°(n)). Furthermore, they satisfy the following properties: 
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• ^ ^ strongly in L^(0, T; L^(ri)) as r —)■ 0. 

• Vi weakly* in L°°{0,T; L°^{k})) as r —0 for i = 1,... ,n. 

• There exists C > 0 such that for all t > 0 and i = 1,... ,n, 

(27) ||C’'’^^||L 2 ( 0 ,'r;/fl(C)) + ^llL 2 ( 0 ,r;_H'l(ll)) +T ^WVi ^ — '^rVi ^\\L^{T,T-,H^{ny) < C, 

where TTrrl[\-,t) = ri^^\-,t — r) for t <t <T is a shift operator. Let D C M"’ be a compact 
domain such that p^'^\x,f) = ... ,pn^){x,t) G D for a.e. {x,t) £ Ll x {0,T). Then, for 

all f G (7^(1);M"’), up to a subsequence, as r —^ 0, 

^W/(^(r)) ^ ^W/(r^W) strongly in L\t),T-L\LI)). 

Since and {rjf'*) are assumed to be bounded in L°°(0,T; L°°{kl)), the strong conver¬ 

gence also holds in LP(0,T;LP(i2)) for all p < oo. This theorem extends [191 Lemma 13], 
proved for /(s) = s, to arbitrary continuous functions /. 

Proof. The proof is based on the compactness result in m Lemma 13], whose proof goes 
back to [3, and an induction and approximation argument. We perform the proof in two 
steps. In the first step / is assumed to be a monomial, in the second step we approximate an 
arbitrary continuous function by a polynomial and apply the Stone-Weierstrass theorem. We 
set Qt = n X (0, T). 

Step 1. Let f{r]) = r/" := • • • 77 “", where a = (ai,..., an) G Nq is a multiindex. The 

proof is an induction argument on the rank |q:| = cii > 0 of the multiindex. If |a| = 0, 
the statement is trivially true. Let us assume that ^ ^ 77 " strongly in L?'{Qt) 

as r —)• 0 for all a G Nq with |a| < fc, A: > 0. Let a G Nq be a multiindex such that 
|a| = A: -|- 1 > 1. Then there exists an index zq G {1,..., tt,} such that ai^ > 1. Hence, we can 
define the multiindex /3 such that fij = aj — for j = 1,..., ti and |/3| = k. 

Introduce = ^N)(r 7 N))^ and y = . Clearly, is bounded in L°°(0,T; L°°(H)). 

Since the multiindex jS has rank k and thus satishes the induction assumption, y 

strongly in L/^^Qt). We claim that and vf^) are bounded in L^(0,T; H^{Tl)). 

Indeed, it follows from (f27|l that is uniformly bounded in 

L'^{Qt)- As a consequence, 

VyM = +.e^"V(77^) 

= (^(N)/3veW+ ^ 

fc:/3fc>0 \i¥^k j 

is uniformly bounded in L?'{Qt), and (y*-^^) is bounded in L^( 0 , T; LA^(H)). In a similar way, 
we can show that (y^’"^ r][J) is bounded in L^(0, T; LA^(H)). Applying [T9l Lemma 13] to the 
sequences (y*-’^^) and {v^), we infer that there exists a subsequence, which is not relabeled, 
such that UVio strongly in L‘^{Qt), which means, by dehnition of and /?, that 

^W( 7 ^ 0)/3 ^fjh strongly in L\Qt). 

Step 2. It follows from the previous step that the statement of the lemma is true if / is a 
multivariate polynomial. Let / G (7^(11; M"') be given. Since D is compact, we may apply the 
Stone-Weierstrass approximation theorem to obtain, for any e > 0, a multivariate polynomial 
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P : D ^ such that \f{ri) — P{r])\ < e for r/ S -D. Since and ^ are bounded in L°°, 

we have for some C > 0 , which does not depend on e, 

W&'fiP) - < Ce, l|{F(x) -{/Wllt^iQrt < Ce. 

Thus, 

\\Pf{P)-a{ri)\\L^(Qr) < W&^fiP) - PPi'iPWiJiQr) 

+ 1\&'P{P) - (PmLHQr) + ll«P(l) - UiMLHQr) 

Since P is a polynomial, the hrst step of the proof applies and the last term on the right-hand 
side converges to zero as r —)• 0 (at least for a subsequence), resulting in 

limsup )-?/(??) IIl 2 (Qt) ^ 2 Ce. 

r —>-0 

Since e > 0 is arbitrary and the left-hand side does not depend on e, it must vanish, finishing 
the proof. □ 

Lemma 8 (Generalized Aubin lemma II). Let be a sequence of functions which are 

piecewise constant in time with constant step size r > 0 and which satisfy a < u^'^\x,t) < b 
for a.e. {x,t) G H x (0,T) for some a, 6 G ffi. Furthermore, let Q G C^([a, 6 ];M”) be a 
nonnegative increasing convex function and assume that there exists C > 0 such that for all 
r > 0, 

- 7r^U^^^||L2(r,T;Hl(D)') < C. 

Then there exists u G L‘^{0,T] H^{Q)) such that, up to a subsequence, 

^ u strongly in L^{0, T; L^{Q)) for all p < oo. 


This result generalizes Theorem 3a in m, stated for Q{s) = s'” with m > 0 . A related 
result has been proved in [22\ Theorem 1]. Instead of the bound on it is assumed 

that the function \Q'\ is bounded from below by a positive value near Too and that the set 
{x : Q'{x) = 0} is finite. Thus, our result seems to be complementary to that one in [ 22 ] . 


Proof. Let (f £ X := P[^ (Q) Cl L°° (Q) be a test function. Then the positive and negative parts 
of 4> satisfy (/>+ = max{0, (p}, 4>- = min{0, 4>} G X. By the convexity of Q, we obtain 

— [ — 7rT-u^'^^)Q'{u^'^^)4>+dx > — [ (Q{u^'^^) — 7TT-Q{u^^^))(j)+dx, 

r Jn -rdn 

— f {u^'^'' — TrT-u^'^'^)TrT-Q'{u^'^'^)4>-dx = — f — u^'^^)Q'{TTTU^'^^){—(j)-)dx 

'T Jq. t Jn 

> - / (TTrQ{u^'^'>) - Q{u^'^'>)){ — (p-)dx = - [ (Q(u‘''"^) - TrrQ{u^'^'>))4>-dx. 

'T Jq t Jn 

Adding both inequalities and taking into account that cp = (p^ + </>-, we find that 
1 


(28) 


T Jn 


— TrrQ{u^^^))(pdx 


1 


< — / — TrTU^'^'l)[Q'{u^'^^)(p-^-+ TrrQ'{u^^'l)(p-)dx 
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< -T:rU^'^'>\\Hi{ny\\Q\u^'''^)4>+ + T^TQ'{u^'"'’)(t>-\\m{n)- 

We estimate: 


^ (|Q'(nW)lV+ + |g'(nM)V0+ + </.+VQ'(n("))|2)dx 

<2 [ |Q'(iiW)|2(4 + |V(/.+ |2)dx + 2 [ (j)l\VQ'iu^^^)\‘^dx 
Jn Jn 

< c(i + ||Q'(K:’':)||;/i(n)) lli^llx- 

In a similar way, we can verify that 




Thus, ([28]) gives 


where 



7TrQ{u^^'>))cj)dx<CF^-\t)U\\x, 


FM(t) = i||(nM -7r.uW)(t)||^i(^),(l + ||Q'(nM(t))||Hi(^ + ||vr.g'(uM)|bi(^,)). 

This means that 

r-i||(Q(nM) -7r.Q(uM))(t)||x' < CF^W. 

The assumptions of the lemma imply that is bounded in L^{t,T). Thus, we obtain 

a uniform estimate for — in L^{t,T]X'). Because of the bound of 

Q{u^'^'>) in L^(0 ,T;and the compact embedding ^ L?‘{Vl), Aubin’s lemma in 

the version of m yields the existence of a subsequence, which is not relabeled, such that, 
as r —)• 0, Q{u^'^^) —>■ Q* strongly in L^(0, T; L^(n)) and a.e. in 11 x (0,T). Since Q is 
strictly increasing, this shows that —> Q~^{Q*) a.e. in 11 x (0,T). We 

set u := Q~^{Q*)- Then the L°° bound for and the a.e. convergence yield —>■ u 

strongly in L^’(0, T; for all p < oo. □ 


3.4. Further results. We show that the Fisher information |Vy^p(l/i is subadditive, and 
we recall a convex Sobolev inequality. 

Lemma 9. Let p be an absolutely continuous measure with respect to the Lebesque measure, 
and let f, g : Q ^ [0, oo) be measurable, bounded, positive functions such that ^/J, yfg e 
H^{Lt,dp). Then 

[ \^VTT^\^dp< [ \Vy^\^dp+ [ |VV5lV 

J Q, J Q, J Qi 

This result was proven in [231 Section 3.6] in a slightly different context. For the convenience 
of the reader, we present the (short) proof. 

Proof. We dehne the function F : [0,1] —)• M by 

F{s)= f \Vy/f\‘^dp+ [ \Vy/sg\‘^dp- [ \Vy/f + sg\‘^dp, 

J ^ J ^ </ Q 


s G [0,1]. 
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Then F{0) =0 and F'{s) > 0 for all s € [0,1] since 

F'{s)= [ [ VyTT^-vf ^ 

Jfi Jfi V 

J ^ t/ Q 


Vf + sg 


dfj, 


-VV/ + sg^ 




V Vf + sg f + sg 

f |Vv/g|^d/i+ [ -jT—\^\/f F s9?dix-2 f ,/f V• VV/+ sffd/i 

Jn Jn J + sg Jq f + sg 


Vf + sg 


dg, > 0. 


We conclude that -F(l) > 0 which shows the lemma. □ 

Lemma 10. Let fl C (d > 1) be a convex domain and let g G he a eonvex function 
such that Ijg” is concave. Then there exists cs > 0 such that for all integrable functions u 
with integrable g{u) and g"{u)\Vu\^, 


1 

W\ 


/ g{u)d: 
Jn 


X- g 


1 

W\ 


udx < 


cs 

|fl| 


f g"{u)\Vu\‘^dx, 

Jn 


where |n| denotes the measure ofLl. 

A proof can be found in [3l Prop. 7.6.1] or [H Remark 3.8]. 


4. Proof of Theorem [T] 

We divide the proof into several steps. 

4.1. Time discretization and regularization of system ([1]). We recall the definition of 
the entropy variable w = h'{u) for u € where h is dehned in ([9]). Lemma [5] shows that h' 
is invertible, thus we may define u = {h')~^{w) for w G R"" and we may set u{w) = u. By 
Lemma [ 6 l the matrix B{w) = A[u){h")~^{u) is positive dehnite for all tc G R and u = u{w). 
We introduce a time discretization for ([T|). Let T > 0, G N, and let r = T/N be the time 
step size. Furthermore, let 0 < e < 1 be a regularization parameter and let m G N be such 
that L°°{Ll) compactly (i.e. choose m > d/2). Given G we wish 

to hnd G L7™'(fl;R"') which solves the discretized and regularized problem 

(29) — f {u{w^) — u{w^~^)) ■ (pdx + f V4> : B{w^)Vw^dx + Vbi;{(p,w^) = 0 

X Jn Jn 

for 4> G L 7 ™(fl;R"’), where 

(30) beiVw^) = [ + ^J)-^w^)dx + e ^ D^cf-D^w^dx, 

^ 2<|o|<m 

and 11“ is a partial derivative of order |a|. We prove the existence of weak solutions to (l29|). 

Lemma 11. Let (jTD-dS]) hold and let '. +L ^ be measurable such that h{u^) G L^(n). 
Then there exists a sequence of solutions G L 7 ™(fl;R"’) to ([ 2 ^ satisfying the discrete 
entropy inequality 

(31) [ h{u{w^))dx + T [ Vw’^ : B{w'^)Vw^dx + T\{w’^,w^) < [ h{u{w’^-^))dx. 

J VI J VI J VI 
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Proof. The idea is to apply the Leray-Schauder fixed-point theorem. Let y G and 

1 ] G [0,1] be given. We first solve the linear problem 


(32) 


a{w, (f) = F{(t)) for all (f G M"), 


where 


a{w,(j)) = / V4> : B{y)Vwdx + T‘^be{w,4>) 
Jn 


F{4>) = -- f {u{y) — u{w^ ^)) ■ (fdx. 


r Jn 


The forms a and F are bounded on The matrix B{y) = A{u{y))h"{u{y)) ^ is 

positive semi-definite, 


v^B{y)v = [h''{u{y)) h”{u{y))A{u{y))[h''{u{y)) ^v] > 0 


for all V G M”, thanks to (I22p . Hence, the bilinear form a is coercive: 


a{w,w) > for w G M"'). 


Therefore, we can apply the Lax-Milgram lemma to infer the existence of a unique solution 


w G to (p2]) . This defines the fixed-point operator S : x 

[0,1] — S{y,r]) = w, where w solves ([32]). 

It holds that S'(y, 0) = 0 for all y G Furthermore, standard arguments show 

that S is continuous (see e.g. the proof of Lemma 5 in [H]). It remains to prove a uniform 
bound for all fixed points S{-,y) in L°°(H;M”). Let w G L°°(H;M"') be such a fixed point. 
Then w solves (f32]l with y replaced by w. With the test function cf = w, we find that 

(33) {u{w) — u{w^~^)) ■ wdx F / Vw : B{w)'Vwdx + T^be{w,w) = 0. 



n 


The convexity of h implies that h{x) — h{y) < h'{u) ■ (x — y) for all x, y £ Choosing 
x = u{w) and y = u{w^~^) and employing h'{u{w)) = w, this gives 



Taking into account the positive semi-definiteness of B{w), we infer from (1331) that 



This yields an Lf"* bound for w uniform in y (but not uniform in e and r). By the Leray- 
Schauder fixed-point theorem , we conclude the existence of a solution w G Lf™(n;R"') to 


(1321) with y replaced by w and y = 1. 


□ 


We derive some a priori estimates uniform in e and r. In the following, we set = u{w^) 


for A: > 1, where {w^) solves dSp. 

Lemma 12. Under the assumptions of Lemma MR there exists a constant C > 0 such that 
for all e, T > 0, 
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fc „ k „ 

+ At Pod + T^'^beiw^ ,w^) < / h{u°)dx, 

j = l ^ J —1 ^ 


where po and 6 are defined in 

Proof. By Lemma [TTl the sequence (w^) satisfies (1291) . Then, taking into account the identity 
: B{w^)Vw^ = Vvf : h"{vf)A{u^)V we deduce that 

f h{u^)dx + T f : h"{u^)A{u^)Vu^dx + T^be{w^,w^) < f h{u^~^)dx. 

J Q, J Vt 

Resolving this recursion yields 

P k ^ k p 

/ h{u^)dx + '^^T / : h"{u^)A{u^)Vu^dx + ,w^) < / h{vP)dx. 


j=i 


i=i 


Then the conclusion follows from Lemma [6] and | ^^1=1 


□ 


4.2. The limit e —)• 0. Let {w^) be a sequence of solutions to (1^ . We fix A: G {1,... ,n} 
and set (f = 1,..., n + 1) and (i = 1,... , n). The identity 

{B{w^)Vw% = {A{u^)Vu% 
shows that solves 

(35) - f {u^ - u^-^) ■(j)dx + '^ j 

'T Jn Jn 

- ■Vfidx + r'^beiw^A) = 0 

for all (j) = (i^i,..., fn) G iL™'(n;M"'). We wish to pass to the limit e —)• 0 in ([35]l . 

By Lemma [12] and definition (1301) of 6e, we have 


(36) 


er^ ^ lIrcLI? 
i=i 


1=1 


where here and in the following, (7 > 0 denotes a generic constant independent of e and r. 
Thus, because of the boundedness of {h")~^ (see Lemma [5]), 

I|Vu(^)||l2 (^) = ||(/i"(u(^)))-iVu;(^)||i2(^) < C||Vu;(^) 11^2(0) < 

Together with the L°° bound for (tt^^^), this implies that 

Therefore, up to subsequences, as e —)• 0, 

^ u weakly in W{^), —)■ u strongly in L‘^{Q) and a.e. in Q, 

since H^{Q) embeddes compactly into L‘^{Q). We infer that u^^+i ~ Sr=i Hn+i := 

1 — Ui strongly in L?‘{PL) and a.e. in 12. The L°° and bounds for (u^^^) as well as the 
bound for in (l34|) show that 
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is uniformly bounded in and hence, 

We employ the a.e. convergence of and and the continuity of pi and q to obtain 

-t UiPi{u)q{un+i)^^'^ a.e. in 0, 

and, by the dominated convergence theorem, strongly in L‘^{Q,). Thus, using the bound, 

u[‘'^Pi{u^^'^)q{u^^l^)^^^ UiPi{u)q{un+i)^^^ weakly in 
Similar arguments, using the uniform estimates coming from (j34h . show that 

(37) q{u'n+i)^^‘^ strongly in L?{Q) and weakly in 

(38) ^ g(ttn+i)^/^rt|'^^ weakly in 77^(17). 

It follows from the bound (l36l) that, up to subsequences, 

—)• 0 strongly in i7”^(fl), w weakly in 

We set := u. The above convergences holds for all /c = 1,..., 77, where T = Nt. Thus, 
we obtain a sequence of limit functions {u^). The above convergence results are sufficient to 
pass to the limit e —)• 0 in ([3^ . resulting in 


(39) 


Jn Jn 

- 3u{pi{u)q{ui^^^)^^'^Vq{ul^^^)^^‘^] ■ V(f)idx + f {w ■ (j) + Vw : V(f))dx = 0 


for (j) G 77™'(n;M”). By density, this relation also holds for all 4> G 77^(11;M”). Note that 
generally we cannot identify w with {h')~^{u) anymore but this is not needed in the remaining 
proof. 

Finally, we wish to pass to the limit e —?■ 0 in ()34l) . where has to be replaced by 
Since 


(40) 


q{u 




the strong convergence (ri|^^)^/^ —)> in L^(H) and the weak convergences (1371) and (l38]l 

imply that 


9(«i+i)^/^V(uJ''V/2 -^V{q{un+iY^'^uJ‘^) - uJ\q{un+if/‘^ 
= g(un+i)^'^^Vu|'^^ weakly in L^(n). 


In fact, since by ((Ml) 


I k (^^Si) K'^ 11L2 (H) 


< Cr-i/2 
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the above weak convergence also holds in L^( 0 ). In particular, by the weak lower semiconti¬ 
nuity of the norm, 


lim inf 
£—^0 


> [ q{un+i)\Vuy'^\'^dx, 
Jn 

liminf [ |> f \Vq{un+i)\‘^dx, 

1 Jn 

liminf II11^1 > ||u;||^i(f^). 


£—^0 


£—^0 

Recall that = u and = w. Passing to the limit inferior e —)• 0 in 
that > llw?*-^^ll^i(Q)) we infer that 

(41) f h{u’^)dx + 4 tpoY^ / 9«+i)^|V(uj)^/^|^dx 

Jn Jn 

k „ k „ 

-|- Arpod E/ |Vg(un-i-i) ^ I dx + T ||'R’'^'||ffi(c2) ^ / h{u )dx. 

j = \ ^ j = l 


and observing 


4.3. The limit r —^ 0. We set u^'^'>{x,t) = u^{x) and w^'^'>{x,t) = w^{x) for a; G H, t G 
{{k — 1)t,/ cr]. Equation (p^ can be formulated as 

-[ [ -7TrU^^^)-(j)dxdt + '^ f j 

Jt Jn ^ Jt Jn 

(42) - 3u(^Vi('w^^^)g(^‘i+i)^'^^Vg(u„+i)^/^] • V(f)idxdt 

+ [ f ■ (p -|- : Vp)dxdt = 0 

J T Jo. 

for all p{t) G if^(n;M"') being piecewise constant in time and, by density, for all p G 
L^(0,T;if^(n)). Inequality (^1]) becomes 

r pT p ^ 

/ h{u^'^\T))dx + 4po / / V |V(ufV/^|^dj;(it 

Jq j 0 Jfl i=\ 

+ 4po6 [ [ \Vq{u'^\Pi^/‘^\^dxdt + T‘^ f ||i/;^’'^||^i(j^xdt < f h{vP)dx. 

Jo Jn Jo ^ Jn 

This gives the following uniform estimates: 

(43) lk(?^i+l)^^^V('u(^V^^IlL2(0,T;L2(D)) + lk(«i+l)^'^^ llL2(0,T;LIl(f^)) - 

(44) ''■||'li’^’'^||L2(0,T;Hl(f^)) - 

(r) 

These bounds as well as the L°° bound for (u) '^) show that 

+ qiu^nhy^'^^2 

j=l ^ 3 / 

= u'ppi{u^^^)Vq{u2lif/'^ 
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q{u 


i'r) 

n+1 




is uniformly bounded in L^( 0 , T; L^(n)) and consequently, 

(45) ll«rVi(w^^^)g(wl+l)^^^llL2(0,T;/fl(O)) < C. 
Similarly, (H3|) yields the estimate 

(46) \\{ufV^W:Um\LH».T ;HHn)) < C. 

Thus, the L°° bound on and estimates (H3]l and give 


-111 (t) ('r)ll 

T Wu'' > -'KrU'^ 


(47) 


< 


^lk(«2l)^^^IU“(r,T;L-(n))||V(ufVi(^^^^^)9(«2l)^^^)IL2(^_r;L2(o)) 

2=1 

n 

+ 3^\\u^^Piiu^^'>)q{un+if^^\\L °°(r,T;L°°(0)) l|Vg(«n+l)'/'|| L2(0,T;L2(a)) 
2=1 


Now, we dehne the function Q{s) = q(ay^‘^da for s G [0,1]. Then Q G C'^([0,1]) is 
nonnegative, convex, and strictly increasing. It holds (see (1431) 1 


(48) 


||(5^(^ti+l)||L2(0,r;_ffi(O)) < C. 


By assumption ([7]), /q 'is uniformly bounded a.e. and thus. 


vg(^l:|i) 


Q"(«Si) 


9'(wi+i) 


is uniformly bounded in L^(0, T; L^(n)). We conclude that 


(49) IIQ(«i+l)llL2(0,T;Hl(O)) <C- 

Estimates (I47l) - (l49p show that the assumptions of Lemma [8] are fulfilled, and we infer the 
existence of a subsequence, which is not relabeled, such that, as r —)• 0 , 

(50) ul^^i ^ Un+i strongly in L'’(0, T; L'’(n)), r < oo. 

This result, the bound ([43l) . and the continuity of q imply that 

(51) g(un+i)^/^ strongly in L'’( 0 ,r;L'’(fl)), r < oo, 

(52) ^ qiun+i^^'^ weakly in L^(0, T; 27^(11)). 

Using the L°° bound for (uj’^^), we have, up to a subsequence, Ui weakly* in 

L°°(0, T; L°°(n)) as r —)• 0. This convergence also holds in L^. Thus, (l50l) implies that the 
relation = 1 — is satisfied by the limit function, Un+i = 1 — X]”=i Ui. The set 

{v G L2(0,r;L2(L!)) : u > 0 a.e. in 11 X (0,T)} is (strongly) closed and convex. Hence, it is 
also weakly closed, and the property > 0 holds in the limit, i.e. Uj > 0 a.e. in 11 x (0,T). 
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It) 

We turn to the convergence properties of the sequences {u^ ) for i = 1,... ,n. We cannot 
expect strong convergence of but the generalized Aubin-Lions Lemma [3 shows that the 

product converges strongly, where / is any continuous function. To make 

this precise, we verify the assumptions of LemmaEl Set := and := u\'^\ 

Because of the L°° bounds for up to a subsequence, 

Vi = weakly* in L°°{0,T;L°°{Q)). 

Furthermore, by (1^ . ^ C = Q{un+i)^^'^ strongly in L^(0, T; L^(i7)). Estimates (l43]l . 

ii, and (I47p show that the assumptions of Lemma [3 are satisfied, and we conclude the 
existence of a subsequence (not relabeled) such that 

fiv)^ = f{ui)q{un+i)^^‘^ strongly in L^{Q,T;L^{9.)) 

for any function / G We choose f{s) = and f{s) = SiPi{s) for s = (sj) G 

Then 

{ui^f^‘^q{u'nlif^‘^ u]^'^q{un+if/‘^ strongly in ^^( 0 , T; L 2 ( 0 )), 

(53) UiPi{u)q{un+if/‘^ strongly in L^(0, T; L^(0)). 

We conclude from the bounds (j45p and (j46p that the above sequences converge weakly in 
L^(0,T; and the limit functions can be identified: 

(54) ^ qiun+if^"^ weakly in L^(0, T; iL^(fi)), 

(55) ^ UiPi{u)q{un+if^’^ weakly in L^(0, T; iL^(n)). 

We infer from estimate (I47p that 

T~^{ui^ — TT-rul'^'^) dtUi Weakly in L^(0, T; L7^(n)'), i = 1,... ,n. 

Moreover, taking into account (I44p . 

—)• 0 strongly in L^( 0 , T; 77^(0)). 

These convergence results as well as the convergences (I51l) - (l53p and (|55p allow us to perform 
the limit r —)• 0 in (1421) . which yields the weak formulation (|15ll . 


4.4. Entropy inequality and positivity. It remains to verify the entropy inequality (I16p 
and the (conditional) positivity of Un+i- Since the entropy density h is convex and continuous, 
it is weakly lower semi-continuous [5l Corollary 3.9]. Thus, by the weak convergence of 




/ 


h(u{t))dx < liminf / h{u^'^\t))dx for a.e. f > 0. 

^^0 Jn 


Employing the convergences (l50l) . ([ST]), and (fMp . it follows that 

converges weakly in but because of the bound (143 p this convergence also holds in 1?“. 

Q'(iii+i)V(ttfV^^ ^ q{un+i)^Ui‘^ weakly in L^(0, T; L^(n)). 

These results, together with (|52p . allow us to pass to the limit inferior r —)• 0 in (I4ip . yielding 
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Finally, assume that 

(56) f \logq{s)\d 

Jo 


s = +00 for all 0 < 6 < 1 . 


We deduce from the discrete entropy inequality (HD) and definition Q of h that 

f f logq{s)dsdx < f h{u^'^\x,t))dx < f h{u^)dx for a.e. t > 0. 

Jq Ja Jo. Jq, 

Then, by the strong convergence of and the nonnegativity of /^logg(s)ds > 0, 

we can apply Fatou’s lemma yielding 


If 


U-nJf-l (lT,£) 


logq{s)dsdx < / h{u^)dx. 

Jn 


In particular, log ( 7 (s)ds < oo for a.e. x € Q. We conclude from this fact and 

assumption ([56|) that Un+i{x,t) > 0 for a.e. x € 0 and t € (0,r), which ends the proof. 

5. Proof of Theorem [3] 

We define the relative entropy density 

(57) h*{u\u°°) = h{u)-h{u°°)-h\u°°)-{u-u°°) foruGM^. 

We split h* in several parts, /i* = /i^ + /i 2 + / 13 , each of which is nonnegative, where 


i=l 


hi (u|u°°) = ^ I Ui log -Ui + u. 




h*^{Un+l\u°°) = 


“"+1 q{s) 9(^C+l)„.oo 


run-\-i 


log- 


log- 


u“+idcj. 


q{u^+i) Ji Qiu^+i) 

n 

hl{u\un = x{u) - y(u-) - - ur) logPi(u-), 

i=l 

where y is dehned in ([8|). The entropy inequality (I16p and the conservation of u{t) give 

n n-^ p TL 

(58) / h*{u{t)\u°^)dx + co I / (q{un+if'^\Vu]^‘^\^ + \Vq{un+if^‘^\^Idxds 

Jo. Jo Jn ^ ^ 

< f h*{u^\u°°)dx, t > 0. 

Jn 

We prove now that the above entropy inequality, reduced to an inequality for h^, and the 
convex Sobolev inequality in Lemma [TOl yield exponential convergence of Un-\-i{t), while the 
entropy estimate for hi and the logarithmic Sobolev inequality allows us to conclude the 
convergence of Uj(f) for i = 1 ,..., n. 

Step 1: Exponential convergence of Un+i{t). Let g{s) = log da for s € [0,1]. 

This function is convex since g''{s) = /q{suff^i) > 0 by assumption. Again 

by assumption, 1/g" = q/q' is concave. Choosing (/>j = 1 in the weak formulation 
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(USD and summing the equations from i = it follows that /u^j^^dx = 

= |fl| for t > 0, and in particular, 


9 


mL 


^n+l 


dx = g{l) = 0 . 


In “n+l 

Thus, we may apply the convex Sobolev inequality in the version of Lemma [TOl 


J^h*2{un+i\u°°)d: 


“n+l , 

<X = -r^ / g 


Jn \u^+iJ 

CS f q'{Un+l), 


f-n^^\dx<^-^ f g" 


/ Un+l \ 

Jn"" Ku^+J 


V 


Un+l 


-''n+l 


dx 


Jn 


Jn q{Un+l) 

By assumption, q' is strictly positive on [0,1], i.e. 0 < gi < q'{s) for s G [0,1], so 

CS f q'iun+i)"^ 


Therefore, 


yields 


^ / h*2{Un+l\u°°)dx < f J-^JjJ^^\VUn+l\^dx 

fI Jn 9i|“l Jn q{un+i) 


qi\^\ Jn q{un+i) 

4cg 

gi|Q| 


\Vq{un+i)^/^\^dx. 


J. 


h*2{Un+l{t)\uX^^)dx + 


cpqi 

Acs 


h 2 {un+i{t)\u’^_^i)dxds < I h*{u^\u°°)dx, 


0 Jn 


and Gronwall’s lemma gives 

(59) [ i)(ix < [ h*{u°\u°°)dx. 

Jn Jn 

The strict positivity of q' implies that the function s i—)• is strictly convex. More¬ 
over, = 0 and = 0. Therefore, by a Taylor expansion, 

/i 2 (ttn-i-i|'U^i) > (7/2)(tin+i — Inserting this inequality in (l5^ gives (fTTp . 

Step 2: Convergence for {ui{t)). We assume that q{s) > go > 0 for s G [0,1]. It follows 
from the entropy inequality (|58p that 


f hl{u{t)\u°°)dx + coqo f f ^ |< f h*{u^\u°°)dx, t > 0. 

Jn Jo Jne Jn 

We apply the logarithmic Sobolev inequality on bounded domains with constant cl > 0 m 
Lemma 1], 

« n « n « 

/ h\{u{t)\u°^)dx = '^ j Ui log JJ^dx < CL \Vu]^‘^\^dx. 

Jne 7i Jci 7 ^Jn 

& 2=1 ‘ 2=1 

Inserting this inequality into the entropy estimate gives 

f hl{u{t)\u°°)dx +[ hl{u{t)\u°°)dx < f h*{u^\u°°)dx, t > 0, 

Jn Cl Jne Jn 

and then, Gronwall’s lemma shows that 


/ hl{u{t)\u°°)dx < / h*{u°\u°^)dx, t > 0 . 

In Jn 
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Finally, since = |u°°)| = 0, and d‘^h\/duiduj = 5ij/ui > 5ij for u G S’, 

we obtain h'[(u\u°°) > |u — which proves estimate (1181) and hnishes the proof. 


6. Proof of Theorem [U 

Let u = (til, • • •) Un) and t; = (t>i,..., Vn) be two bounded weak solutions to ([I])-([ 2 ]). Since 
Pi = 1 for all i = 1 ,..., n by assumption, dU becomes 

(60) dtUi = div {q{un+i)Vui - UiVq{un+i)), i = 

Summing these equations from t = 1,..., n, the equation for Un+i = 1 — Yl7=i reads as 

(61) dtUn+l = div {q{Un+l)VUn+l + (1 - tin+l)Vq'(tin+l)) = AQ{Un+l), 

where Q{s) = fQ{q{cr) + (1 — a)q'{a))da for 0 < s < 1. Furthermore, VQ{un+i) • p = 0 on 511, 
t > 0 and tin+i(0) = ti°_|_i := 1 — similar equations holds for Vn+i- Since Q is 

a nondecreasing function, we can apply first the method to (1611) to show uniqueness for 
the (n + l)th component, i.e. Un+i = t>„+i. Second, we employ the convexity of the entropy 
to prove that tij = u* for i = 1 ,..., n. 

Step 1: Uniqueness for tin+i- Let t > 0 and let C{t) £ H^{Q) be the unique solution to 


— AC{t) = {un+i—Vn+i){t) in n, VC-P = 0 on H. 

We know that Un+i — Vn+i G L^(0, T; L^(H)). Thus, t i—>• C{t) is Bochner integrable and 
C G L^{0,T;H^{Q)). As dt{un+i — Vn+i) G L^(0,T;Ff^(H)'), we have even the regularity 
AdtC G L^(0,r;iL^(r2)'). Therefore, using (f 6 T]) . we obtain for a.e. t > 0, 


In 1^'’^ {-^dtC, 0 = {dt{un+i - W+i), C) 

= - / V{Q{Un+l) - Q{Vn+l)) -VCdx 
JQ 

— / i^Qi^n+l) (5(r’n+l)) (rtn+l Vn+l)d: 

JQ 


Here, (•,•) again denotes the duality pairing of iL^(H)' and Ff^(H). The right-hand side is 
nonpositive since Q is nondecreasing. This implies that 


[ |VC(t)pdx < [ |VC(0)|2dx, t > 0. 

Jn Jq 

At time t = 0, —AC(0) = {un+i — u„+i)(0) = 0 in H, thus VC(0) = 0. Hence, |VC(t)| = 0 a.e. 
in n, which gives {un+i — Vn+i){t) = —A^it) = 0 in H. 

Step 2: Uniqueness for (tti,... ,Un)- Let 0 < e < 1. Similarly as in [16], we introduce the 
distance 


deiu.v) (^f.eiui) + f,eivi) - 2 ^^ 

where ^^(s) = (s + e)(log(s + e) — 1 ) + 1 , s > 0 . 


dx, 


As fs is convex, we have Cei^i) + feivi) — 2^e((uj + Vi)/2) > 0 in H and hence, de{ui,Vi) > 0. 
We need the regularization e > 0 since Ui and Vi are only nonnegative and thus, expressions 
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like log((ttj + Vi)/2) may be undefined. Since Un+i = Vn+i by Step 1, we may abbreviate 
q := q{un+i) = q{vn+i)- Then, using (jHO]) . we compute 


d ^ ( 

—4(u, = X] ( + e)) + {dt'Vii log(uj + e)) 

2=1 ^ 


- ( dtiui + Ui),log 


Ui + Vi 


+ s 


^ / f (gVui - UiVq) ■ + {qVvi - ViVq) ■ 

^ In V Ui + e Vi + e 


- {qV{ui + Vi) - {ui + Vi)Vq) 


Viuj + Vj) 

Ui “t“ T 2s 


dx. 


Rearranging the terms, we arrive at 


n « 

•si. 

n p 

+E 


Vui\'^ \Vvi\'^ _ |V(ut + 


lij + e Vi s Uj + Uj + 2e 

Ui + Vi 


i=l 


Uj _ 

Ui+ s Ui + Vi + 2e 

Vi Ui + Vi 


n \<Ji 


Vi + e Ui + Vi + 2e 


Vq ■ Vuidx 
Vq ■ Vvidx 


n « 

- 4 ^ / {\Vy/ui + e\‘^ + \Vy/vi + s\^ - \V\^Ui + Vi + 2 ep) qdx 

i=i 

+ 2 ^/ (—-f- ^ y/qVy/q-Vuidx 

+ 2^/ (—-f- y/qVy/q-Vvidx. 

'^,jQ.\Vi + e Ui + Vi + 2eJ 


2=1 


Now, we apply Lemma [9] with djjL = qdx and f = Ui + e, g = Vi + e, showing that the first 
integral on the right-hand side is nonnegative. We observe that ^^(^(O), u(0)) = 0 as u and v 
have the same initial data. Thus, integrating the above expression in time, we obtain 


Ui -I- Vi 


Ui £ Ui Vi 2s 


(62) deiu{t),v{t)) <2'^ f f ( 

Jo Jn V 

^ ft f 

-h 2 V / / , 

Jo Jn \'^i ^ Ui + Vi + 2s J 

Since Vyig, yJqVui, y/qVvi € L^(0, T; and 

Ui + Vi 


y/qVy/q ■ Vuidx 


\ y/qVy/q-Vvidx. 


Ui + s 


< 1 , 


Vi+£ 


< 1, 


Ui -\- Vi -\- 2s 


< 1 , 
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the dominated convergence implies that the right-hand side of (1621) tends to zero as e ^ 0. 
From the nonnegativity of we deduce that d^{u{t),v{t)) —)• 0 as e —)• 0, which means that 

( U ■ “ 1 “ V ■ \ 

" ^ M 0 ase-S^O a.e. in x (0, oo). 

According to Taylor’s formula, there are functions 6^, % : n x (0, oo) such that 


= e 


Uj + Vj Uj - Vj 




+ fC' { Oe- 


Ui + Vi 


+ (1 


^eiVi) = Ce 


Ui + Vi Ui- Vi 


= ie 


Ui + Vi\ (Ui + Vi\ Ui- Vi 1 „ ( Ui +Vi 


-ie 


+ + (1 





Adding these identities and employing the estimate C”{s) = {s + s) ^ > 1/2 for 0 < s < 1, 
we infer that 

^eiui) + - Vif. 

This estimate and ([631) prove that Ui = Vi in Q x (0, oo) for i = 1,..., n. 


7. Extensions 


In this section, we discuss some extensions of the diffusion system ([T]). 
Reaction terms. Cross-diffusion systems with reaction terms. 


(64) dtu — div(A(u)Vu) = f{u) in 11, t > 0, 

can be treated similarly as in m- More precisely, if there is a constant cj > 0 such that 
f{u) ■ h'{u) < Cf{l + h{u)) for all u £ then there exists a global weak solution to ([2]) 
and (|64p . The proof proceeds as for Theorem [H where the right-hand side of the entropy 
inequality (|3H) has to be replaced by 

/ h{u^)dx + T / f{u)-h'{u)dx< / h{u^)dx + rcf / {l + h{u))dx. 

Jn Jn Jn Jq 

Then, for sufficiently small r > 0, the integral rcf h{u)dx can be absorbed by the left-hand 
side of (j34p . For instance, reaction terms of Lotka-Volterra type 


fi{u) 



i = 1,... ,n, Sij> 0, 


are admissible. The large-time behavior result is valid only under an additional condition on 
/(n), namely f{u) ■ h'{u) < 0 for tt G If we suppose conservation of the “total mass”, i.e. 
Sr=i = 0) the method allows us to prove uniqueness for Un+i. Uniqueness for the 
remaining components Ui follows if there exists C > 0 such that for all Ui and Vi, 


E 




‘2'iii 

Ui + Vi 


+ fi{v) log 


2vi 

Ui + Vi 


2 = 1 




Ui + Vi 


2 
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where ^(s) = s(logs — 1) -|- 1 (see the proof of Theorem S]) . More general conditions on f{u) 
can be found in m- 

Drift terms. In the presence of environmental or electric potentials or of chemotactic signal 
concentrations, the diffusion system contains additional drift terms, 


(65) 


dtu — div(^(u)Vu + D{u)V4>) = 0 in If, f > 0, 


where D{u) = {Dij(u)) is an n x n matrix and the ith component of D{u)'V4> is given by 
Dij{u)V4>j, where 4>j = 4>j{x) is some potential. Assume that h is such that Vu : 
h”{u)A{u)Vu > gi{u)\V for some nonnegative functions gi{u). Then, using the test 
function h'{u) in the weak formulation of (j65p . we compute 


dt 




Vu : h”{u)A{u)Vudx 



Vtt : h”{u)D{u)V(f)dx 



n „ \ ^ f 

/ gi{u)\S/ui\^dx + - Y / 


Gk{u)\V(l)k?dx, 


where we employed the Cauchy-Schwarz inequality and have set Gk{u) = Y17j=i9ii'^) 
Djk{u)‘^ with H = h"{u). Thus, if Vcfi is bounded in and Gk{u) in L°°, we achieve 
some gradient estimates, which are the basis for the existence analysis. An example is the 
ion-transport model [8] 

Aij{u)=Ui forz/j, Aii{u) = Ui + Un+l, Dij{u) = UiUn+Aj 


for z, j = 1,..., n. The entropy density can be dehned by 

n 

Hu) = Y {uiifogUi - 1) +Ui4>i) + Un+l(logUn+l - 1). 
1 = 1 


Then the Hessian h"{u) does not depend on <f>i. A formal computation, using the Cauchy- 
Schwarz inequality and the identity YHi=i ^Ui = —Vun+i, gives 


fL 

dt 



h{u)dx 





V 



Uj 

Ufi+l 



2 

dx 


< - 


E 

i=l ■ 


/ 1 


V log 


Ui 




- |V0i|2 dx 


! L n I O n 

Y / [‘^Un+l\yuY\^ + \S/Un+l\'^ +2\VuYi?)dx+ Y / «i«n+l|V(/)i 


2 = 1 


2 = 1 


As q{s) = s in this model, we hnd the same estimates as in the proof of Theorem [T] (also see 
[71 Section 3.2]). This shows that our strategy can be adapted to cross-diffusion systems with 
drift. 

Other diffusion coefficients. Our main assumption on the transition rates is that they are 
given by the product ofpi(u) and qi{un+i) (see Appendix!^. Also other choices are possible. 
An example is the diffusion system of [24] . which is derived from a stochastic lattice model 
by assuming that the transition rates are given by Pi{u) + qi{un+i) for some special functions 
Pi and qi. The diffusion matrix has the structure 

A(„\ - fo‘iH-U2) + U2 (ai-l)ui \ 

{a2-l)u2 02(1 - ui) + uiJ ' 
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where ai, > 0. The corresponding diffusion system possesses the entropy density 
2 

h(tt) = '^Ui{logUi - 1) + (1 - Ul - M 2 )(log(l -Ui - U2) - 1), u= {ui,U2) G 


2=1 


and the new diffusion matrix B = h"{u) ^A{u), given by 

/(ai(l - ui - tt2) + -U 1 U 2 

I -U1U2 {a2{l - ui - U2) + Ui)u2) ' 


B = 


is symmetric and positive semi-definite on For our analysis, we need bounds from 
h”{u)A{u) (see LemmaE]), which are less obvious since 


Vulh''{u)A{u)Vu 2 = 


ai|(l - ri2)Vui-I-niVu2p a2|tt2Vtti-|-(1 - tti)Vu2| 


ui{l - Ul- U2) 
4|V/uIu^|^, 

only yielding an bound for V-^uiU 2 in L^. 


+ 


^2(1 - Ul - U2) 


Appendix A. Formal derivation of the n-sPECiES population model 


We derive formally the cross-diffusion system ([T]) from a master equation for a discrete- 
space random walk in the diffusion limit. We consider random walks on a one-dimensional 
lattice only, since the derivation can be extended in a straightforward manner to the higher¬ 
dimensional situation. The lattice is given by cells Xj (j G Z) with the uniform cell distance 
h = Xj — Xj-i > 0. The proportions of the zth population in the jth cell at time f > 0 is 
denoted by Ui{xj) = Ui{xj,t). The species move from the jth cell into the neighboring cells 
j ± 1 with the transition rates T/’ . The master equations are given by 

dtUi{xj) = T^~^'~^Ui{xj-i) -h T^~^^~Ui{xj+i) - (r/’+ Tj~)ui{xj), i = l,...,n, 
and the transition rates are defined as 

n 

(66) T/’ = aoPi{u{xj))qi{un+iixj)), Un+i{xj) = I-'^Uk{xj), 

k=l 

where u = {ui,... ,Un)- The quantities Pi{u{xj)) and qi{un+i{xj±i)) measure the tendency 
of the species i to leave the jth cell or to move into the jth cell from one of the neighboring 
cells, respectively. More precisely, Ui{xj) denotes a volume fraction of occupancy and Un+i 
the volume fraction not occupied by the species. Our assumption is that the transition rates, 
measuring the occupancy and the non-occupancy, separate, resulting in the product of pt and 
qi- Other choices are possible (see [23] for an example), but the analytical treatment of the 
corresponding diffusion systems is not obvious. 

For the derivation of the diffusion model, it is convenient to introduce the following abbre¬ 
viations: 


P- = Pi{ui{xj),... ,Un{xj)), qj = qi{Un+l{Xj)), 

9kPi = ^{uiiXj),...,Un{Xj)), dqj = qi{Un+l{Xj)). 
OUk 


Thus, we can rewrite the master equation as 
(67) 


fJo ^dtui = ql {pI ^ + pI^^uI^^) - pIuI (g^+^ + ql ^) 
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Set D = dx- We compute the Taylor expansions of pi and qi {i = 1,... ,n) and replace 
— u^, by the Taylor expansion zizhDu^j^ + ^h?D^ui + 0{h^). Then, collecting all terms up 
to order 0{h?‘), we arrive at 


J di 1 j 

Vi =Vi 




+ hY,dkPlDul + ^ dkplD‘^u{ + dlfplDulDul + 0 {h? 


k=l 


,k=l 


k,l=l 




Qi =qI^ hdplDul^^ + — {dqlD^Y + ^^qI (^<+ i )^) + Oih^) 


= qj^ hdqj Y + 
k=l 




-dqj Y^^^i + d^Qi Y I +0{h^). 

k,i=l 


k=l 


In the last step, we have used Un+i = 1 — Yl^=i '^k- We insert these expressions into (|67)) and 
rearrange the terms. It turns out that the terms of order 0(1) and 0{h) cancel, and we end 
up with 


Uq ^dtul=Y (qM ^ik + qI ul dkPl + Pi ul dql) 

k=\ 

n 

+ {2qldkpl6ii + qluldl^pj - pjujd'^q’). 


k,l=l 

We choose cjo = h~‘^ and pass to the limit h —>• 0: 


” ( dPi \ 

dtUi = \ qiPi5ik + Qi'^i^:^ +PiUiq'ij 

f dPi 9^Pi 

+ Y DukDue + Qi'^i a„. al. - Pi'^iQi 


k,e=i 


' dukdul 


A lenghty but straightforward computation shows that the last sum equals 

n 

DukD 

k=l 

and we end up with 

n 

dtUi = dY, Duk 

k=l 

which is the one-dimensional version of ©• 


QiPi^ik T Qi'^i 




dpi 

duk 


■PiUiQi 



QiPAk + QiUi^^ + pmq'i 
OUk 
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